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Using the mixed space representation {t,p ) in the context of scalar field theories, we prove in a 
simple manner that the Feynman graphs at finite temperature are related to the corresponding zero 
temperature diagrams through a simple thermal operator, both in the imaginary time as well as in 
the real time formalisms. This result is generalized to the case when there is a nontrivial chemical 
potential present. Several interesting properties of the thermal operator are also discussed. 
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I. INTRODUCTION 



In a series of recent papers Q, |M |^, it was shown for theories involving scalar as well as fermion fields that 
every graph in momentum space at finite temperature in the imaginary time formalism ^ |^ is related to the 
corresponding graph of the zero temperature Euclidean field theory through a thermal operator which has a rather 



' simple form. Namely, for a scalar iV-point amplitude (at any loop) at temperature T one has 



e)m=n(l + n,(l-S.)). (2) 

i=l 

Here / characterizes the total number of internal propagators and V denotes the total number vertices in the graph 
(with the usual relation for the number of loops L = I — V -\-\), rii = n{Ei) corresponds to the thermal distribution 
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associated with the internal propagator carrying energy Ei = \Jk^'^ + m? and Si — S(Ei) is a reflection operator that 

, changes Ei — s- —Ei (namely, it gives a term with Ei —Ei). We denote the internal and the external three momenta 

^ ' of a graph generically by {k,p ) respectively and 5''v'\k,p) enforces the appropriate three momentum conservation at 
the vertex v. For simplicity, we have included in a delta function which reflects the overall conservation of the 

external three momenta. Furthermore, 7]^ represents the zero temperature graph in momentum space obtained 
' from the Euclidean field theory. This remarkable result is, of course, calculationally quite useful since the worrisome 
sum over the internal discrete energy values (particularly at higher loops) has already been reduced to evaluating zero 
temperature energy integrals. More than that, this allows us to study many questions of interest at finite temperature 
such as Ward identities, analyticity 0, B| more directly. The original proof of this result in momentum space 0, Q 
however, is quite involved and uses regularization procedures that obscure the origin of such a relation. Furthermore, 
the proof leaves one with the feeling that such a relation is particular to the imaginary time formalism. In thispaper, 
we discuss a simpler proof of the thermal operator representation both in the imaginary time formalism [j, |5|, |g| as 
well as in the real time formalisms 0]. Furthermore, we also extend this relation to the case when there is a nontrivial 
chemical potential ^ and point out various interesting aspects of this relation. 

At finite temperature, it is already noted that simplifications arise when one works not in the energy-momentum 
space, but rather in a mixed space where energy has been Fourier transformed 0,0,^3- exploit this feature to 
give a simpler derivation of the thermal operator representation in both imaginary time and the real time formalisms 
with and without a chemical potential. In this paper, we will discuss in detail theories involving scalar fields only 
to bring out the essential underlying features. The remarkable feature of the thermal operator relation, in such a 
representation, is that while the finite temperature graph depends on 



iP ^lTiT.E,,T^), (3) 

where Ta,a = 1,2,-- - ,N represent the external time coordinates of the graph, the thermal operator depends only 
on (T,Ei) but not on the external time coordinates (the zero temperature graph depends on (Ei,Ta) but not on 
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T). As a result, the time derivative operator (with respect to external times) commutes with the thermal operator 
and the discussion of our results holds equally well for theories with fermions as well as Yang-Mills fields which we 
will discuss in detail in a separate publication. Our paper is organized as follows. In section II, we first discuss 
the thermal operator representation in the imaginary time formalism without a chemical potential and then with a 
nontrivial chemical potential. In this section, we also point out various properties of the thermal operator which are 
quite interesting. In section III, we derive the thermal operator representation in the closed time path formalism 
[a, [Til Il2l | , again with and without a chemical potential, where the general proof is really much simpler than the 
imaginary time formalism. In section IV, we discuss the thermal operator representation for a general time contour, 
that includes the case of thermofield dynamics where the thermal operator, in general, is not as simple as in 

the imaginary time and the closed time path formalisms. We conclude with a brief summary of our results in section 
V. 



II. IMAGINARY TIME FORMALISM 



Let us consider a massive real scalar field theory in Euclidean space. In this case, we know that the zero temperature 
propagator in momentum space is given by 

A^^-'Hpe,E) = = (4) 

where we have defined 



The energy variable can now be Fourier transformed to give the propagator in the mixed space at zero temperature 
to be 

A(^=o)(r,i?)^ I^^Z^^J-[eir)e-^^ + 9{-r)e^^], -oo < r < oo. (6) 

At finite temperature in the imaginary time formalism (we will set the Boltzmann constant to unity for simplicity), 
the propagator in the momentum space has the same form as in Q with pe ^ 2TTkT where k is an integer. In this 
case, the Fourier transform of the propagator leads to 

A(^)(r, E) ^ ^^^g^ ^ ± [{0{r) + n{E)) e'^^ + (0(-r) + n{E)) e^^] 

= ^ [^W { (1 + niE))e-^- + niE)e^^} + e{-T) {n{E)e-^^ + (1 + n{E))e^^]] , -1 < r < 1, (7) 

which is symmetric under t —t — r. It is important to recognize that, in the imaginary time formalism, time is 
rotated to the negative imaginary axis and lies between the interval (0, ^) so that the propagator (being a difference 
of time coordinates) is defined only within the interval (— y, y). Furthermore, let us note that, in this mixed space 
representation, the thermal propagator Q can be naturally written as the sum of a zero temperature part and a finite 
temperature part much like in the real time formalisms 6] 

Am(r,S) = A(^=o)(r,S)+A(T,i;), (8) 
where the finite temperature part has the form 

A(r,i?) = ^(e-^^ + e^O- (9) 

Furthermore, we note that both the zero temperature propagator and the finite temperature propagator Q satisfy 
the same equation 

(^^-E^^A{T.E) = -5{r), (10) 
but the finite temperature propagator satisfies the periodicity condition 

A(^)(r<0,i?) = A(^)(T+^,i?), (11) 
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following from the KMS condition It follows, therefore, that A in ((Sj) satisfies the homogeneous equation 



i?^ A(t,£;)=0, 



(12) 



and is responsible for incorporating the periodicity condition Hll|l . 

It follows now from the forms of the propagators in ^ and (jT)) that we can write 



- O(^)(i?)A(^=0)(T,i?), 



1 1 

— < T < — , 

T - - T' 



(13) 



where we note that the basic thermal operator 0'''^\E) is independent of the time coordinate of the propagator. 
As we will see, this basic factorization of the finite temperature propagator in terms of a thermal operator (that 
contains all the temperature dependence, but no time) and the zero temperature propagator (which carries all the 
time dependence) is at the heart of the thermal operator representation of any finite temperature graph. Given the 
factorization in (|13|l . it immediately follows that any one loop graph with N external lines (see figure ^) would lead to 
the thermal operator representation (we consider the (f)^ theory for simplicity, neglect the overall coupling constants, 
assume that all momenta are incoming and identify kN+i = ki,pN+i = Pi) 




FIG. 1: One loop diagram in the 4>'^ theory with A'^ external time coordinates. 
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=1 

where we have identified 

N N 



FIG. 2: Two loop diagrams in the (f)^ theory without internal time coordinates. 



Similarly, the thermal operator representation immediately follows from the factorization of the propagator H13|l for 
any higher loop graph where all the vertices have only external times (no internal vertices present). It is obvious, for 
example, in the case of the graphs shown in figure [51 in the (j)^ theory. 

The difficulty in establishing the thermal operator representation for an arbitrary graph arises when there are 
internal vertices present for which the internal time coordinates have to be integrated over all allowed values. For 
example, in the 0^ theory the self-energy graph at two loops can have diagrams containing internal vertices of the 
form shown in figure |3| (r, f in this case correspond to internal time coordinates that have to be integrated over the 
allowed values). Keeping in mind our earlier comments, we note that at finite temperature, the time integration goes 





FIG. 3: Two loop self-energy diagrams in the theory with internal time coordinates r, f . 



over 



dr, 



while in zero temperature graphs, the internal time needs to be integrated over 



dr, 



(16) 



(17) 



and it is not clear a priori how the range of integration IjlGI) can be extended to H17|l in order to establish the thermal 
operator relation. Furthermore, the finite temperature propagator is only defined in the interval (— i,i), and it 
is not clear if it can be analytically continued to other regions (preliminary analysis indicates it cannot be done so 
consistently). Before giving a general proof that such an extension of the range of integration can, in fact, be made, 
let us work out a simple example to bring out some of the essential features. 

Let us first consider the product of N propagators with a common time that is being integrated, namely, 



N 



(18) 



where r^, z = 1, 2, • • • , iV are assumed to be external times that lie between the interval (0, y). Using H13|) . we can 
write the above expression also as 



N 



4=1 

N 



N 



h = X{{l + n,{l-Si)) / dr[]A(^=o)(r-r„ii:, 



z=l 
N 



/oo ^* 
dT\[A(^='>\T-n,E,)--h, 
-oo -1 



(19) 



5 



where we have defined 

N 



1=1 



(20) 



We note here that it is because the basic thermal operator in the factorization of the propagator is independent 
of the time coordinate that it can be taken out of the time integraL Furthermore, although the finite temperature 
propagator cannot be extended beyond its domain, once we have extracted the thermal operators, the zero temperature 
propagators are defined on the entire real axis, a fact which we have used in the above. Since the external times 
satisfy < < using the definition of the zero temperature propagator in © the integrals in the bracket above 
can be evaluated in a simple manner and lead to 





Tii e^'^' 



N 



n 



-EiTi 



■i=i ' / ^« * 1=1 * 
It is now straightforward to show that the thermal operator annihilates H21|l , namely, 

-E^Ti 



(21) 



^ / ^ \ 1 ^ -E0 

na +».(!- 5.)) i-nTT;^(-^.) E^nv 

i=i \ i=i ' ^ / ^« * i=i * 



(22) 



This, therefore, establishes that for a product of N propagators integrated over a single common time which is 
integrated, we can extend the range of integration to write 



N 



N 



N 



(23) 



(We want to emphasize here that even though the basic thermal operator is independent of time, it is improper to 
take the products of the thermal operators inside the integral in (|23() and write it as a product of thermal operators 
being integrated over the interval (— oo, oo) since the thermal propagators are not defined outside of their domain. 
Such an attempt would lead to various divergences as well as inconsistency problems.) The case when there is a single 
internal time integration works out in a simple manner because the extra terms in (|2(JII do not involve any nontrivial 
time ordering as is clear from H21() . This is no longer the case when there are two or more internal time integrations. 
We can always extend one of the time integrations to the entire real axis as discussed above. But, once this is done, 
the subsequent integrals will involve nontrivial time ordering. Nonetheless, case by case, one finds explicitly that the 
integration range can be extended to the entire real axis (after factoring out the thermal operator) when two internal 
time integrations are involved. This, in turn, suggests that there must be a general proof for such an equivalence for 
an arbitrary number of internal time integrations which we discuss next. 



A. General Proof 



We have already seen that when there is one internal time coordinate that is being integrated, the range of the 
integration can be extended from (0, ^) to (— cxd,c») under the action of the thermal operator. The generalization 
of this result to an arbitrary number of internal times that are being integrated can be carried out in the following 
way. First, let us note that since the zero temperature (as well as the finite temperature) propagator satisfies H10|l . it 
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follows that 

N 



k / J oc 

.00 ~, 
1=1 

TV 

= n'A(^=°)(r.-r„£;,). (24) 



i=l 



Here rfe is any external time coordinate and the prime on the product implies the absence of the term with i — k. As 
a result of this identity, we can write 

pJr\{l^iT^'\r-n^E.)^[~^^+El) ' {\{ A^^='Hr, ~ E..)^ . (25) 

In general, a homogeneous term (a term annihilated by the differential operator) is allowed on the right hand side 
of (|25|l . However, since a Feynman diagram is a time ordered quantity, a homogeneous term is not expected on 
physical grounds. That this is true mathematically can also be seen as follows. Let us evaluate the integral (|25() in 
the momentum space. Using the definition of the Euclidean propagator in momentum space in we obtain 

/ drnA(^=°)(r-r„i?.) 



27r („2 +^2^ 



i=l 



The delta function allows us to eliminate one of the piE variables and let us choose the dependent variable to be pkE 
related to r^. Eliminating this variable, we obtain 



/ drnA(^=°)(r-r„i?0 



N N ■ , s 

^ ^ I TT' ^P^E TT' e-^^^-^--^^ 

-1 JV 



which is the result obtained in We would hke to note here that the operator ^—^3 + -^l) can be thought of 



7 



in terms of the standard integral representation as 



92 



El 



da 



E 

1 



■,-aEi 



Nl 



-a 



.N 



1 + E 



N=l 



1 

Md4 



N' 



(28) 



The relation H25|) is quite useful in proving that the range of the finite temperature integration can be extended to the 
entire real axis under the action of the thermal operator for any number of internal time integrations. The important 
thing to note here is that the differential operator commutes with the thermal operator (since it is independent of 
time) and that the equivalence can be established recursively as follows. First, we note from H23|l that for a single 
time integration this is true. Let us assume that we have a product of propagators with two internal times that are 
integrated. The most general form for such a product can be written as 



N2 



N, 



drdf n^^'^^(^-^-^'') n^^''^(^-^'^") \\[^'-^\r-f„E, 



\i=i 



0(T) 
0(T) 



N2 



drdf n^^''^''^(----^^) \\{^^^='\r-f,E^)\ \\[/^^^='\f-f,^E,) 



No 



N3 



dr 1^ df(l[M^="HT^n,E,)j (^nA(^=o)(r-f,^„)j j^J] A(^=o) (f - f„ i?,) 



92 



\j=i 

-1/ Nl 



N2 



O^^^ (~^ + E!) in A(^-o) (n ^r..E.)]l df ( n A(^-°) {r, - f, E.)][ U A(^-°) (f - f„ , E, 



N3 

n 



where iVi, N2, N3 are arbitrary integers. Here we have identified 



O 



(T) 



^Ni \ / / 

n (1 + n,(i - s,)) n (1 + ""(1 - '^") 



N3 



(29) 



(30) 



as well as used H25|l with t\ as the external time coordinate for simplicity. This shows that the case of two internal 
time integrations can be reduced to that of a single time integration where we know from H23|) that the range of 
integration can be extended to the entire real axis (under the action of the thermal operator) so that we have 



'JVi 



JV2 



A'3 



drdf (nA(^)(r-r„i?,)) ( J] ^^''H^ " j (^JI A(^)(f - f„, S,) 



0(T) 
0(T=O) 



i=l 

dr- 



\a=l 



'1/ Nl 



J + [n (-1 - ^^)] fj^ [ n A^^=°^ in-f, E^)j n a^^=°^ ir-f^,E,) 

j drdf (^nA(^=°)(r~r„i?,)j (^1] A(^=°) (r - f , ^„) j j^J] A(^=°) (f - i?,) j , (; 



where we have used ()23|l to restore the t integration. This process can be used recursively to show that the range 
of integration can be extended to the entire real axis (under the action of the thermal operator) for any number of 
internal time integrations. This, therefore, proves the thermal operator representation for any arbitrary graph with 
N external legs, namely. 



/ n H n (^^fe\KP)iT - / ri H n (2-)^^^^^^,^)^ 

•' i=l ^ ' v=l '' i=\ ^ ' 11 = 1 



(T) (T=0) 

!N 



(32) 
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There are several things to note here. The proof of the thermal operator representation in the mixed space is more 
direct and the origin of this relation can be traced to the factorization of the thermal propagator in terms of the 
basic thermal operator which is independent of time and the zero temperature propagator. There is no necessity for 
classifying the graphs into trees or for introducing any regularization as one does in the momentum space analysis 
(which is quite unusual since the finite temperature results are not expected to be divergent). As we will show later, 
the derivation of the thermal operator representation in the mixed space is even simpler in the closed time path 
formalism. We also note here that although this analysis seems to suggest that this equivalence holds only for graphs 
with external legs, such a relation holds even for graphs without any external legs which follows in a straight forward 
manner from the closed time path formalism to be discussed later. Here we simply note that if we are looking at the 
pressure in the (p^ theory at two loops at finite temperature (see figure ^ we have the explicit result (neglecting the 
overall factors involving the coupling as well as the symmetry factor) 



El , kj 




FIG. 4: Pressure diagram in the (j)^ theory at two loops. 



where 

Jo 

J_ 

= O(^) dTdfA^^="\T-f,E,)A(^=°\T-f,E2)M^="\T-T,E3) 

Jo 

The thermal operator for this graph is given by 

3 

= J] (1+71,(1-50) • 



(33) 



(34) 



(35) 



As we have discussed earlier, in the case of a single integration over time, the range of integration can be extended to 
the entire real axis under the action of the thermal operator. Using this as well as shifting variables of integration, 
we obtain 



jiT) = o(T) 



{2Ei)i2E2)i2E3) ^ ,12 3; 



oo Jo 



/' df / dr^i^i±|^i::^A(^=°)(r,i?i+i?2 + i?3) 



{2Ei)i2E2){2E3) 

^ OmLs,,_^o r dr 2(^1 + E2 + E,) e,+E. + E,), 

T " ' J-oo {2Ei){2E2){2E3) ^ ' " 

where we have used the basic definition of Fourier transform for a finite interval (fco is an integer) 



(36) 
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Recalling that in the continuum limit 

^4o=o,o ^ 2^5(0), (38) 

we immediately identify that 

J(T) ^0(T)j(T=O)^ (39) 

and the thermal operator representation works even for graphs without any external legs. This is more directly seen in 
the closed time path formalism that we will discuss later. It is worth noting here that in the imaginary time formalism, 
the graphs without any external leg always have a factor ^6ko=o,o which must be identified with the continuum case 
as in H38I) . This is not necessary in the real time formalism where time is a continuous variable defined over the entire 
real axis. 

Let us also comment here on some interesting aspects of calculations in the mixed space. The calculations in 
this case are more like the real time calculations in the sense that the amplitudes contain all possible factors of the 
statistical distribution function. However, when the results are Fourier transformed into energy-momentum space, 
because of various identities, the number of statistical distribution functions in an amplitude reduce to one per internal 
loop 1','^. We note that given the thermal operator representation (|32|) of a graph in the mixed space, one can go to 
the energy-momentum space in the standard manner. Here, we have to simply remember that the zero temperature 
amplitude j^^~^\Ei, r^) is a function of external times which are restricted to lie between < < i. As a result, 
the Fourier transform needed to go to the energy-momentum space is that over a finite interval (involving integer 
energy) even though we have a zero temperature amplitude in the Euclidean space. 

The thermal operator representation for any graph at finite temperature is a remarkable result. Physically, a 
Feynman graph at finite temperature represents an ensemble average while a zero temperature graph corresponds 
to a vacuum expectation value. Therefore, a relation between the two can exist only if the expectation value for a 
string of operators in any complete set of states (say in the energy eigenbasis) would be proportional to the vacuum 
expectation value of the same string of operators. Although at first sight this seems unlikely, let us show that this is 
plausible with the simple example of the propagator for a massive, real scalar field at the tree level (which will also 
explain the factorization for the thermal propagator) . Using the standard field decomposition in the Euclidean space, 
we note that we can write (in the mixed space) 

0(t, P ) = ^ (e-^^ a{p ) + at {-p )) , (40) 



where as before, E = y^p'^ + m?- At zero temperature, this leads to 

(O|r(0(ri,pi)0(r2,P2)) |0) = (5(3)(p- + _i_ [e{T^ - ^^y-E,(r^-r.) ^ ^(^^ _ (41) 

On the other hand, in any eigenstate of energy containing N quanta of momentum p, we have the expectation value 

(7V,p|T(0(ri,pi )0(r2,P2 ))\N,p) 

= _1— (jV,p I [0(ri - T2) {a{n)a\-p2)e-^-^^+^-^^ + a\-n)a{p2)e'^'''-''^'') 

+e(T2 - Ti) {a{p2)a\~n)e''^-^-^--- + a\-p2)a{p^)e-^^^^+^-^-)] \N,p) 

= 5(3) (pi +p2)5(3)(p'_p^) J_ \e{T^ _ T2) f(l + iV)e-^^(^i-^^) + Ne^'^-^'-^A 

ZEi L V / 

+9{t2 - Ti) (^Are-Ei(^i-^2) + (1 + 7V)e-^i(^i-^^)) 

= <5(3)(pi +p^)S(^)(f-p^) (1 + iV(i _ S{Ei)) kn - r2)e-^i("i-"=) + 0(r2 - ri)e^i("^-"^)" 

2Ei L 

= <5(3)(p~pi) (l-fiV(l - SiEi)) (O|r(0(Ti,pi)^(T2,p2)) |0). (42) 

This shows that the expectation value of the time ordered product of two fields in any higher energy state is propor- 
tional to the expectation value of the same operators in the vacuum state and the proportionality factor is reminiscent 
of the basic thermal operator. In fact, carrying out the thermal ensemble average, this proportionality factor indeed 
becomes the basic thermal operator of ()13|l. 
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B. Some Properties of the Thermal Operator 

The basic thermal operator that leads to factorization of the thermal propagator has several interesting features. 
In this section, we discuss some of them that are relevant to a better understanding of this factorization. This will 
also be quite useful in connection with the study of factorization in the real time formalisms. Let us note some of the 
basic properties of the reflection operator. By definition {S{E))'^ — 1 and for the bosonic distribution functions we 
have (this is identical to the well known result n{—E) = —(1 + n{E)) for a bosonic distribution function) 

S{E)n{E) = S{E)^— = ^^S{E) ^ -{l + n{E))S{E), 

S{E){1 + n{E)) = S{E)-l^ ^ -^-L^siE) = -n{E)S{E). (43) 

gT — 1 e T — 1 

Using these basic properties, it is easy to see that the basic thermal operator can be written in various ways as 

0^^\E) = 1 + n{E){l - S{E)) = (1 + S{E)){1 + n{E)) = -(1 + S{E))n{E)S{E). (44) 
Furthermore, it follows from the definition that 

{o'^'^\E)^^ = {l + n-nS){l + n-nS) 

= (1 + nf - n{l + n) + (n^ - n{l + n)) S 

= {l + n-nS)=0'^^\E). (45) 

The basic thermal operator, therefore, is a projection operator. Consequently, the inverse of this operator does not 
exist and the thermal operator representation for graphs does not have an inverse relation. We note from H44I) that 

S{E)0^^\e)^0^^\e), (46) 

so that we have 



1 - S{E) 



2 y > ^ 2 

Let us note that we can define another operator 



(47) 



0^^\E) = l + n{E){l + S{E)), (48) 
which also corresponds to a projection operator and satisfies 

io'^\E)f = o'^\e), ^±tmio^-\E) = 0, ^±^^0'^\e) = 0'^\e). (49) 
The two projection operators, however, are not orthogonal and satisfy 

0^^\e)O^^He) = (1 + 2n{E))O^^HE), 0^^\e)o''^\e) = (1 + 2n{E))0^'^\E) . (50) 

The thermal operators, of course, depend on the temperature. Denoting the temperature dependence explicitly, we 
can write 

0^^\E)^l + n{E,T){\^ S{E)). (51) 

It can now be directly checked that 

C)(T2)(£;)0(Ti)(£;) ^ 0^^'\E). (52) 

At first sight this may seem a bit strange. However, this result is quite consistent with the underlying physics. Let us 
recall that the effect of the thermal operator is to reproduce an ensemble average. Once the ensemble average has been 
done through the operator O'^-^^-' (i?), the resulting amplitude is a scalar (proportional to the identity operator). The 
application of a second thermal operator O^'^^^E) to the result is then equivalent to thermal averaging the identity 
operator which simply results in a multiplicative factor of unity. 
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Finally, to understand the meaning of the thermal operator O^^^ (E) as a projection operator, let us note from (jSj) 
as well as (|13|) that 

O^^He)A^^=°Ht,E) = A(^)(r,i?) = A(^=")(r,i?) + A(T,i?). (53) 

We know that the zero temperature propagator does not satisfy the periodicity condition (|ll|l . Rather, it is the 
temperature dependent term A(t, E) that enforces the periodicity condition. Thus, we can think of the thermal 
operator as projecting on to the space of functions satisfying the periodicity condition. Of course, it follows from the 
definition that 

0'^^\e)A{t, E) ^ 0, (54) 

which is consistent with the fact that A(r, E) is the homogeneous solution of the Green's function equation H1(J|) 
(recall that the thermal operator commutes with the differential operator). The meaning of the projection operator 

O (E) is also now clear since it can be directly checked that 

O^^^ {E)A{t, E) = 2(1 + n{E))K{T, E). (55) 



C. Chemical Potential 



So far our discussion has been within the context of a canonical ensemble where there is no chemical potential. 
Chemical potentials arise when we have a conserved charge and we are dealing with agrand canonical ensemble. In 
this case, the Hamiltonian in the definition of the partition function is generalized to |5j 

H ^ H -fiQ, (56) 

where Q represents the conserved charge and /x the chemical potential associated with it. In the context of a scalar field 
we can introduce a chemical potential if we are dealing with a complex scalar field where there is a natural definition 
of a conserved charge. The free Lagrangian density, for such a theory in the presence of a chemical potential, can be 
written as 

C^{dt+ 0* [dt -i^J-)(j)- V0* • V0 - m'^(t)*(t). (57) 

It can be checked that the Lagrangian density (|57|l leads to the Hamiltonian H56|l where Q represents the conserved 
charge that generates the global phase transformations of the system. We note that the addition of a chemical 
potential can be viewed as introducing a constant electrostatic potential into the system. We also recall here that for 
a relativistic massive boson, the chemical potential has to satisfy 

/i < m. (58) 
In Euclidean space, the Lagrangian density (|57|l takes the form 

Ce = {dr + m) 0* {dr - pi)(f> + ^(f>* ■^(p + m^(j)*cf>, (59) 
which leads to the zero temperature propagator in momentum space to be 

where as before we have 

E^=p^ + m^. (61) 
This can be Fourier transformed to give the mixed space representation of the zero temperature propagator to be 



A(T=o^f^){T,E) = J 



dpi 



277 {pE - il^)'^ + E^ 

2E 



^ e{T)e~^^-f'^^ +0i-T)e^^+^^^ . (62) 
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At finite temperature, the momentum space propagator continues to be given by (|60|l with pE = ^irkT where k is an 
integer. Fourier transforming this, we obtain the mixed space representation for the thermal propagator to be 



-2i7TkTT 



2TTkT - ifx)^ + 



= 1^ + n{E - Ai))e-(^-^)^ + (0(-t) + n{E + /i))e(^+^)"' 

ZE L 



1 

2E 



6I(t) {(1 + n_)e-(^-^)^ + n+e(^+^)^} + 6l(-r) jn^e-^-^-^)^ + (1 + n+)e(^+^)^} 
where, for simplicity of notation we have defined 

n± — n(E ± ^). 



(63) 



(64) 



We note here that when ^ = 0, this reduces to {Tj) as it should. However, unlike the real scalar field, here the propagator 
carries a direction, namely, the direction of the charge flow (from (j) to (j)*). The propagator is not symmetric under 
T ^ — T. This is a reflection of the fact that the chemical potential inherently distinguishes between particles and 
anti-particles. However, under the simultaneous reflection (r, ^) — (t, /x), the propagator is invariant. Furthermore, 
both the zero temperature and the finite temperature propagators satisfy the equation 



{{dr-^,r-E')A^^''^Hr,E)^-S{r) 



(65) 



The finite temperature propagator, however, satisfies the periodicity condition (|ll|l . 

The form of the thermal propagator in the presence of a chemical potential is rather complicated and it is not clear 
whether there will be a factorization in this case. A little bit of analysis, however, shows that the propagator can, in 
fact, be factorized as 



A(^.-)(r,£;) = ^ 



e(r) |(1 + n_)e-(^-'^)" + n+e^^+^'>^} + 0(-r) jn^e-^^"^)" + (1 + n+)eS^+^'>^} 
a(^-^)(i?,9,)A(^="'^)(T,i?), 



where 



1 + ^^^^(1 - S{E)) + ^^+^(1 + S{E))^ (a. - f^) 



(66) 



(67) 



We note that the first two terms are similar to the ones in the thermal operator of the earlier section while the last 
group of terms in the above relation is new and vanishes when fi = 0. It is this term that reflects the asymmetry 
in the chemical potential for particles and anti-particles. This basic thermal operator reduces to the one in H13I) and 
continues to be independent of the time coordinate. However, we would like to point out a further simplification that 
takes place in this case. 

Let us note from H62|l and H63() that the dependence of the chemical potential in the exponents of the propagators 
completely factorizes. 



A(^=o.f^\T,E) = e^^A^^=^=^\T,E), 

A(^''^)(t, E) - e^^^ [e{T) {(1 + n_)e-^^ + n+e^^} + 6{-t) {n^e"^^ + (1 + n+)e^^}] 



(68) 



where A(^=o=^)(r, E) is given in ©. Since e^"^ factors out in both zero as well as finite temperature propagators, we 
can write a simpler factorization for the thermal propagator as 



A(^^^)(r,£;) = e'^^C'(^''^'(£^,90A(^=°='^)(^,£^), 



where 



0'^^^^'\E,^r) = 



1 



^^il-SiE)) + ^ 



^il + SiE))^dr 



(69) 



(70) 



Furthermore, let us note that in any IPI graph involving closed loops, the overall factor e^"^ would cancel and hence 
can be ignored. (There are various ways of seeing this. Since the time comes back to itself in a closed loop, this 
factor reduces to identity. In terms of electrostatic analogy, if the particle comes back to the starting point, there 
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is no change in voltage. Such a simpHfication, however, would not take place in a tree level graph which is not 
IPX.) As we have argued earlier, the thermal operator representation of a graph is a reflection of the factorization 
of the thermal propagator and so given the factorization in (|66|l it would seem that we can write a simple thermal 
operator representation for any finite temperature graph even in the presence of a chemical potential. In general, 
however, because of the time derivative terms in the basic thermal operator, this factorization is not as simple as in 
the case without a chemical potential. In the case of graphs, where every propagator is connected to an external time 
coordinate, the thermal operator representation takes a simple form for any IPI graph involving closed loops, 

/ ri H n eHkM^^^^ - / n H n s^j^Hk,p)o^^-h^r"^^\ (td 

where 

Qi^'t^) = Yl {E, ,±drj, (72) 

i=l 

where Tq is the external time coordinate to which the propagator with energy Ei is connected (± in the time derivative 
represents the phase that may arise in changing the argument to the external time coordinate). This is almost like 
the case when there is no chemical potential. 

However, if there are propagators in a diagram which are not connected to an external time, it is not clear a priori 
whether a thermal operator representation can be written for such a graph. The difficulty arises because if both the 
time coordinates associated with a propagator are internal times, it would seem that the basic thermal operator 17U|) 
for a propagator cannot be taken out of the time integration and, consequently, it is not clear whether a thermal 
operator representation of the graph can result. That such a representation, be it nontrivial, may arise can be seen 
from a simple nontrivial graph like the one shown in figure El 




FIG. 5: A three-loop vertex correction diagram in the (0*0)^ theory with two internal time coordinates. 

Let us consider the vertex correction graph at three loops in the complex 0** theory. The graph involves two internal 
times T, f that need to be integrated over and there are two propagators with time coordinates that are completely 
internal. We consider the graph with the charge flows as shown in the figure. In this case, we note that we can write 



(T,m) 
74 



J drdf A(^^^)(ri - t,E,)A^^^''\t2 - T,E2)A^^'^\f - n, ^i)A(^'^) (f - ra, ^2) 
X a(^'^) (t - f , E)A'■^■^''^ (r - f , E) 

= n0(^'^)(i?„a|f))0(^^'')(^„-a|f)) I drdf A(^=o=^)(Ti-r,i?i)A(^=o=^)(r2-r,£;2) 

,^(T=o=p)(~ _ ^^,^^)A(T=o=M)(^ _ r2,^2)A(^-^)(r - t,E)A^^''^\t -t,E), (73) 



1=1 



where the superscript on is to specify on which propagator the time derivative is going to act. To see how the 
time derivatives in the last two propagators can be taken out of the integral, let us note that we can write 

A(^''')(T-f,S) = (^X{E) + Y{E)^^^A^^=''=^'\T~T,E), (74) 

where we have identified for simplicity 

X{E) = l+'^±±^{l~S{E)), Y{E) = '^±-^{l + S{E)). (75) 
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Using this as well as using the identities 



2(£: + £;) 
{2E){2E) 



A(^=o=t^)^r-f,E + E) 



ia,A(^=o=^)(r - f , E) l9,A(^=o=^)(r - f , 
E E 



E + E 



.E)A 



(T-f,£;) 



we can write 



^(t.m) _ £;) aC^^^) (r - f , E) 

= X{E)X{E) + Y{E)Y{E)\ A(^=°='') (t - f , £;) A^^^o^^) (r - f , i?) 

^ hi -\- hi ^ 



+ 



r,E) 



(76) 



(77) 



The T derivative can now be integrated by parts inside the integral and put on the other propagators where the 
argument of the derivative can be changed to an external time coordinate and can be taken outside the integral. This 
shows that although naively we will not expect a factorization of the graph in figure where there are propagators 
containing only internal time coordinates, a thermal operator representation does exist. However, it is not as simple as 
in the case without a chemical potential and furthermore a closed form of the thermal operator can only be determined 
graph by graph. In view of the above analysis, it seems plausible that such a non-trivial factorization may also hold 
for general diagrams. 



III. CLOSED TIME PATH FORMALISM 

As we have emphasized several times earlier^ the thermal operator representation follows much more directly in the 
real time formalism of closed time path [y, [l^, [131 • Let us recall that in the closed time path formalism, the theory is 
defined in Minkowski space where time is a continuous real variable defined over (—00,00), unlike in the imaginary 
time formalism. Of course, the price one has to pay is to double the number of degrees of freedom, for every field in 
our theory (we denote the real scalar field of our theory by 0+) we add another field of the same kind (We refer 
the readers to 6] for details.) As a result, the propagator acquires a 2 x 2 matrix structure and in the momentum 
space has the form 



A^^lip) A^^ip) 



where, for a massive real scalar field, we have 



Af^(p) - hm l^^ + 27rni\po\)S{p^-m^), 

Af2(p) - 2n{e{-po)+n{\p,\))S{p^~m^), 
aL7(p) = 27:{9{po)+n{\po\))S{p'-m^), 



A'^^ip) = hm ~- '-^— + 2TTn{\po\)6{p'-m'), (79) 

with n(|po|) denoting the bosonic distribution function. The components at zero temperature follow from this to be 

aV^=°^(p) = lim 



A^+r'^ip) = 27r9{~po)Sip^~m^), 
aL+=°^(p) = 27:e{po)S{p^ -m^), 
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These are Minkowski space propagators and the ie prescription in the diagonal elements specifies the choice of the 
contour in the complex energy plane (the two diagonal elements simply correspond to time ordered and anti-time 
ordered propagators). Fourier transforming the energy variable, 



A(i,p) = A(i,i?) = j ^ 



dpo _ 



ipot 



A(p), 



(81) 



where, as before, E = yT^^ + m^, for the zero temperature components we obtain 



Af-)(t,i^) 


= lim 


1 

2E 




Afr"'(t,i?) 


1 

2E 








1 

2E 




A^Ir°\t,E) 


— lim 


1 

2E 





(82) 



Similarly, the Fourier transform of the finite temperature propagator yields the components to be 

1 





lim 


T-{t^E) = 


1 

2E 


^-lit,E) = 


1 

2E 


'^lit.E) = 


lim 

, n 



ie~'^* + (l + n(^))e*^*] , 
n{E))e-'^' +n{E)e'^'] , 



iEt 



iEt 



(83) 



We have carefully kept the it terms in the exponent resulting from the Feynman prescription which are essential for 
the convergence of factors in any calculation. 

Looking at the components of the propagators in the mixed space in (I82II and H83(l we see that there is natural 
factorization of the components of the finite temperature propagator so that we can write 



Af_^(i,i?) Af2(t,i?) 



A 



(T=0) 

++ 



(t,E) A'^^-"\t,E) 



A'_l{t,E) A'_l{t,E) 



(1 + n{E)il - S{E))L{e)) A^^=°^ {t, E) = O'^' {E)A^^=°^ (t, E). 



(84) 



Here L(e) is an operator that takes the limit e — s- in the expression on which it acts. If there is no e dependence 
in the expression, the effect of this operator is that of the identity operator. Thus, we see that there is a very 
simple factorization of the thermal operator in the closed time path formalism where each component of the matrix 
propagator factorizes by the same factor which does not depend on time and is reminiscent of (we recall that there 
is no e dependence in the imaginary time formalism). This is, however, not the case in other real time descriptions 
such as thermofield dynamics as we will discuss in the next section. 



A. General Proof 



Given the simple factorization l|84|l of the finite temperature matrix propagators, the thermal operator representation 
for any A^-point graph now follows immediately. Let us note that any graph in the closed time path formalism is simply 
a product of vertices (both "-|-"and "—"types) and the corresponding propagators with integrations over internal time 
coordinates. Since each component of the propagator has the same simple factorization (|84|l . then it follows that 
for any A'^-point graph (with only "+" vertices or "— " vertices or mixed vertices) where there are no internal time 
coordinates. 
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where 7]^ represents the value of the corresponding Minkowski space graph after the energy integrations have been 
carried out and 

/ I 

O^^^ ^l[0^^\E,) + n,{l - S,)L{e)) . (86) 

1=1 i=l 

This is exactly the same as in the imaginary time formalism. When there are internal time coordinates that need to 
be integrated over, however, the proof of the thermal operator representation in the closed time path formalism is 
much simpler. In fact, note that here time is a real variable defined over the entire real axis independent of whether 
we are at zero temperature or at finite temperature. This is the difference from the imaginary time formalism. As 
a result, the thermal operator representation 185|l continues to hold even when there are internal time coordinates 
that need to be integrated over. (Namely, in this case, we do not have to extend the range of integration and thereby 
avoid the complicated proof of equivalence in extending the range of integration as is needed in the imaginary time.) 
We find this proof of the thermal operator representation by far the simplest. Furthermore, this result holds for any 
N point amplitude including the case when iV = 0. (Namely, we have made no assumption about a graph having an 
external leg.) Therefore, the thermal operator representation clearly holds even for graphs without any external legs 
(such as pressure). We have not been able to show this directly in the imaginary time formalism although we have 
argued (based on our result of the closed time path formalism) and shown in simple examples that it must be true. 

It is worth pointing out some features of the action of the thermal operator at this point. First of all, let us consider a 
multi-loop graph with all external vertices only of "-|-"type. At finite temperature, we know that there can be internal 
vertices of "—"type. At zero temperature, however, the "—"type internal vertices give vanishing contribution jl6| . 
There are several ways of seeing this. The simplest reason is probably the most physical, namely, at zero temperature 
any amplitude is simply given by the original theory (one does not need a doubling at zero temperature). Thus, 
applying the thermal operator representation to such a graph would imply 

7f^'+... = o("'7f+t.^ m 

Since the zero temperature graphs would not involve any "—"type intermediate vertices, this would imply that the 
finite temperature result for such an amplitude can be obtained only from graphs involving "-|-"type vertices. This is, 
however, in contradiction with the known fact that we need a doubling of the degrees of freedom at finite temperature. 
The resolution of this puzzle is interesting, which also clarifies the action of the thermal operator in the following 
way. Namely, even when a zero temperature graph vanishes, one should not set it to zero before applying the thermal 
operator to it. Only after the thermal operator has been applied can the relevant terms be set to zero. This is also 
how new channels of reaction 15] at finite temperature can be seen to arise in the thermal operator representation. 
Let us illustrate this with the two loop example in the cj)'^ field theory shown in figure |^ At zero temperature, this 




FIG. 6: A two loop self-energy diagram in the theory with an internal "— " vertex, 
graph leads to (we factor out overall factors involving coupling and symmetry factors) 

/oo 
-00 
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Since both Ei , E2 are positive, this is clearly zero and is consistent with our observation that at zero temperature, 
graphs with internal "—"type vertices do not contribute. At finite temperature, however, this graph does contribute 
and has the value 

J — 00 

2ni(l + ni)(l + 2713) X 2'k6{Ei - E2). (89) 



' —00 

1 



{2Ei){2E2){2E3) 

If we naively set the zero temperature graph to zero, then clearly there is a contradiction. On the other hand, if we 
apply the thermal operator to the zero temperature result, we obtain 

(1 + ni(l - (1 + n2(l - 52)) (1 + n3(l - ^3)) ^-^) ^2E,){2E2){2E^) ^^^^^^ ^ 

(90) 

If we now use the fact that both Ei , E2 are positive to set the first group of terms to zero and use the properties of 
the second delta function, we obtain the finite temperature result in H89|l . This simple example is quite illustrative in 
understanding the action of the thermal operator. 



B. Chemical Potential 



The simplicity of the closed time path formalism continues to hold even in the presence of a chemical potential as 
we will show without going into too many details. Let us recall that the free Lagrangian density for a complex scalar 
field with a chemical potential has the form given in H57() 



(91) 



In the closed time path formalism, we have to double the degrees of freedom and so introducing the doubled degrees 
of freedom (0_, ((>*_) (we label the original fields as {(p+, (t>+))i we note that the propagator acquires a 2 x 2 matrix 
structure (and carries a direction from to </>*) 



with 



A('^'^)(i,£;) = 



A^_^f\t,E) A^^l^\t,E) 
A^^f \t,E) AL^I^'(t,£;) 



A^+fit,E) = 



A'-^l^\t,E) 



A'L'f\t,E) 



A'-^l^\t,E) 



lim 

e^O 2E 



lim 

e^O 2E 

1 

2E 

'2E 
1 

2E 

Ye 



[n_e-*^* + (l + n+)e^^*] , 
(l + n_)e-'(-^-^)*+n+e*(-^+^) 
[(l + n„)e-'^* + n+e'^*] , 



lim 

£^0 2E 

lim 

£^0 2E 



(92) 



(93) 



Here n± are the distribution fimctions introduced earlier in (|64() . 
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Once again, it is obvious that component by component, we have a factorization of the thermal operator as 



{t,E) ^ e'^''0'^^'^'\E,^t)A 



(T=0=m) 



where i,j — ± and 



{l^S{E))Lie) 



-(1 



SiE))Lie)-dt 



(94) 



(95) 



This simple factorization of every component of the thermal propagator as well as the fact that in the closed time path 
formalism, the range of integration over internal time coordinates is the same at zero as well as at finite temperatures 
immediately leads to the thermal operator representation for any IPI graph. However, as discussed in the case of the 
imaginary time formalism, the thermal operator representation in the presence of a chemical potential is nontrivial 
and more involved because of the presence of the time derivative terms in 0'''^'^\E,dt)- 



IV. GENERAL REAL TIME CONTOUR 



There are two commonly used real time descriptions of finite temperature field theory. In the last section, we have 
already discussed one of them, namely, the closed time path formalism |^ 0, ^| which is quite useful in calculating 
various quantities both in thermal equilibrium as well as out of thermal equilibrium. The other commonly used real 
time formalism goes under the name of thermofield dynamics 0, ll,l| which is inherently an operator formalism and 
is quite useful in understanding various operatorial issues such as the nature of the thermal vacuum and the thermal 
Hilbert space etc for equilibrium systems. It can also be used for calculations in thermal equilibrium although its 
main power lies in understanding operatorial issues. Both these formalisms correspond to specific time paths in the 
complex t-plane. In general, one can define a thermal field theory with a time contour in the complex i-plane of 
the form shown in figure [7| jl^l- Here ct is a constant and has the value in the range < ct < 1. When (7 = the 

Im t 



Re t 



r 



T 



FIG. 7: General time contour in the complex t-plane. 

description of the thermal field theory coincides with the closed time path formalism while for ct = i, the description 
corresponds to thermofield dynamics. In this section, we will show that a simple factorization of the thermal operator 
and, therefore, of a thermal graph takes place only for the cases when ct = 0, 1. In the case of closed time path 
corresponding ct = 0, we have already seen this. Here we will study the behavior of the thermal propagator for a 
general time path in the complex f-plane. 

Let us consider a theory describing a real scalar field. The doubling at finite temperature simply corresponds to 
introducing the fields ((/)i, 02 ) on the two paths labelled 1, 2 (the original field that we start out with is considered to 
be 01 ). For a general time contour as shown in Fig. 4, it can be shown that 



A(^)(p) 



a1P(p) Ar2^(p) 



Ag\p) Ag\p) 



(96) 
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where 



Am 



AT) 



12 ip) - e^AV^^b), 



A(?(P) - Am 



(97) 



(T) 

where A^^^ {p),i,j — ± correspond to the propagators in the closed time path formahsm that we have aheady discussed 
in the earher section. Fourier transforming these components of the propagator in the energy variable, we obtain the 
mixed space propagator for the general contour to be 



AT) 



AT) 

^12 



A, 



AT) 



AL7(t 



We note from ^ that in the hmit T -> 0, 



j;,E) = 


1 




sinh 


E 
T 


^,E) = 


1 




sinh 


E 
T 



sinh (1 — cr) 



sinh 



aE 
T 



E 
T 

Af2 



AT-{t,E) 



{t,E) 



sinh 



aE 
'Y 

sinh( (l-a)| 



AL7(t,i?) 



AL7(i,i?) 



(98) 



AT) 

^12 



AZ\t,E) 



it,E) 



At^'\t,E) 

(T=0) 



A^-i)§ A^^-'\t,E) 



,5,,o AV-^^(i,i?) + <5.4 A^^^"\t,E) 



- A 



(T=0) 



{t,E) + e- 



A^^-°\t,E)+6,^o A'-+"'{t,E) 



A 

(T=0), 



(T=0) 



(i,£;) 



(99) 



This shows that if cr 7^ 0, 1, the off-diagonal components of the propagator vanish (at T = 0) leading to a decoupling 
of the two fields at zero temperature. When a — 0,1, the off-diagonal elements of the components do not vanish at 
zero temperature, nonetheless there is decoupling of the two degrees of freedom. 

Since each component of the propagator in the closed time formalism has a simple factorization given by H84|) , it 
follows from 1)98(1 that for a general time contour we can write 



A(^) [t, E) = {1 + niE)il - SiE))L{e)) A{t, E), 



(100) 



where 



A,^{t,E) - Af-"' 
~Ai2[t,E) ^ 

~A2i[t,E) 



sinh ^ 



sinh (1 



A22(t,E) = A 



sinh ^ 

(T=0) 



sinh 



aE 



A 



(T=0) 



sinh 



sinh (1 — tj) 



aE 
'Y 

E 
T 



(T=0) 



^iT=0) 



(101) 



It is clear from the above result that in the general case, temperature cannot be completely factored out of the matrix 



in a simple manner unless ct = 0, 1. For the case of thermofield dynamics where a 



it is easy to show that 



A(^Ht,E) = {l + n{E){l^ SiE))Lie)) { _i 



^^^(f)^ A(^=o)(t,i?), (102) 



so that, in this case, the basic thermal operator takes a matrix form. As a result, the thermal operator representation 
for any graph in the formalism of a general contour (where cr 7^ 0, 1) is not so simple as in the imaginary time 
formalism or the closed time path formalism. However, it is worth noting that although in such cases there will be 
no simple factorization at the level of individual graphs, the simple factorization will occur for the complete set of 
graphs associated with a given physical amplitude (which follows from the fact that a physical amplitude is the same 
in any formalism). 
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V. SUMMARY 

In this paper, we have systematically studied the interesting question of thermal operator representation for Feyn- 
man graphs at finite temperature. By working in a mixed space (r,p ) (or (t^p)), we have given a simpler derivation 
of the thermal operator representation in the imaginary time formalism. We have traced the origin of such a simple 
relation to the fact that the thermal propagator, in this space, has a basic factorization where the basic thermal oper- 
ator is independent of time. We have also generalized the thermal operator representation to the case where there is a 
nontrivial chemical potential. In this case, although the thermal propagator also factorizes, the basic thermal operator 
involves a time derivative which leads to a thermal operator representation for any graph that is highly nontrivial. We 
have tried to study various properties of the thermal operator and have shown that it is a projection operator which 
projects functions into the space where the KMS periodicity condition is satisfied. We have also shown that there 
is a simple thermal operator representation in the closed time path formalism. The derivation, in this case, is even 
simpler than that for the imaginary time formalism. For a general time contour (including the one for thermofield 
dynamics), however, the thermal operator representation is not so simple as in the imaginary time formalism and the 
closed time path formalism. 
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